Upper limits on the gravitational mass loss of the Galaxy 
and the LIGO burst searches 
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We discuss the relevance, for the search of gravitational-wave bursts, of upper limits on the total 
mass loss of the Galaxy which come from various astronomical observations. For sub-millisecond 
bursts we obtain limits on the event rate, as a function of the GW amplitude, which are stronger 
than the corresponding upper limits set by LIGO in the S2 run. A detection of a burst rate 
saturating these limits, with the sensitivities of present and near-future runs, would imply that, 
with some improvement on the accuracy of astronomical observations of the Galaxy, as foreseen 
with the GAIA mission, it might be possible to detect gravitational waves indirectly from their 
effect on galactic dynamics. 



I. INTRODUCTION 

One of the main target of the existing gravitational- 
wave (GW) detectors are short bursts, with a duration 
between fractions of a millisecond and a few seconds, 
which could originate from astrophysical events. The 
LIGO collaboration has published the results of the S2 
run m, (see also f^j for the S3 run) which extended the 
search to lower values of the GW amplitude, compared 
to previous searches by resonant-mass detectors |^ H, |^ . 
Presently, LIGO is performing a long data-taking run at 
its target sensitivity which, in case of no positive detec- 
tion, should anyway produce a more significant bound on 
the event rate. 

The result of these searches can be presented as an 
upper limit on the event rate of GW bursts, TZ, as a 
function of the strength of the GW signal. The latter can 
be conveniently characterized in terms of the so-called 
root-sum-square amplitude /irss, defined by 
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The purpose of this note is to point out that a bound on 
7?. as a function of /ij-ss can also be obtained from consid- 
erations of galactic dynamics, and it is in fact quite signif- 
icant compared to the bounds that can be obtained with 
existing and near-future sensitivities of GW detectors, 
especially for bursts of short duration, say r ~ 0.1 ms. 

The bound emerges from the fact that, if there is a 
steady rate of GW bursts, the Galaxy has a correspond- 
ing rate of mass loss into GWs. Since it is difficult to 
imagine that the burst rate today is significantly higher 
than in the past, one must consider the cumulative effect 
of this mass loss over a period comparable to the age of 
the Galaxy, and this can have significant consequences on 
various aspects of galactic dynamics. These issues were 
first addressed many years ago ^, Ji |a 13 > and in [l3| we 
reconsidered them, using the present knowledge of galac- 
tic dynamics, and we showed that much more stringent 
bounds emerge nowadays. Let us summarize the results, 
referring the reader to Ref. [T^ for details and more ref- 
erences. The main observations that allow us to put a 



bound are discussed below. 



II. UPPER LIMITS FROM GALACTIC 
DYNAMICS 

(i) Effect of the mass loss on the radial velocity of stars. 
If the Galaxy is loosing mass, stars become less and less 
bound and acquire radial velocities with respect to the 
Galaxy rest frame. This modifies the radial velocity of 
stars, Vrj inducing a so-called X-term, 



^ ARiiin2l + KR, 



(2) 



where the term AR sin 21 is the standard effect due to the 
differential rotation of the Galaxy {A is Oort's constant, 
R is the distance of the star from the sun, and / the 
galactic longitude), and the effect of mass loss is in the 
second term, where K — —M/M, and M is the mass 
of the Galaxy. From the experimental bound on K^ we 
deduced in 11Q| a bound 



M < C»(3O)M0/yr . 



(3) 



(ii) Mass loss and outward motion of the LSR. Rather 
than looking at the K term, i.e. at the expan- 
sion/contraction of the stars within a few kpc from the 
sun, one can investigate whether the local standard of 
rest (LSR) has an overall outward radial velocity, as sug- 
gested by eq. (01 ■ Here, the most interesting information 
comes from the observation of the 21 cm absorption line 
toward the galactic center [Il| , which shows that the gas 
along the line-of-sight has a mean radial velocity with 
respect to the LSR of —0.23 ±0.06 km/s. The absorbing 
material is probably at 1-2 kpc from the galactic cen- 
ter. A radial expansion due to mass loss predicts a radial 
velocity Vr ^ r, see eq. (|2l, where r is the radial dis- 
tance from the galactic center, and therefore we should 
expect a difference in velocity between us and this gas, 
Awr = (— Af/Af )Ar, where Ar ~ 6 kpc is the distance of 
the Sun from the absorbing gas. 

In general, there can be both positive and negative 
contributions from different physical mechanisms to the 



value Awr = —0.23 ± 0.06 km/s, and to extract a bound 
on mass loss to GW we require that no fine tuning be- 
tween different contributions takes place. We set conven- 
tionally at 20% the maximum fine tuning that we allow, 
which means that we say that a positive contribution 
from GWs to Avr, if it exists at all, must be smaller 
than 0(0.04) km/s. Setting the distance between us and 
the gas to Ar = 6 kpc, this gives a bound 



M < 0(0.5) Mo/yr. 



(4) 



Of course, precise numbers depend on the level of fine 
tuning that we can tolerate, but it is clear that we can- 
not raise this bound by, say, one order of magnitude, 
without invoking very unnatural cancellations between 
completely unrelated phenomena. 

(iii) Upper limits from globular clusters. Similar 
bounds have been found using globular clusters as 
probes |3,[i2|- The idea is that, if the mass of the Galaxy 
was much bigger in the past, the orbits of globular clus- 
ters would have been much closer to the galactic nucleus, 
and this close interaction with a very massive central nu- 
cleus would have produced the tidal disruption of the 
cluster. The analysis of five different globular clusters 
gives the result 
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M < 0(10) Mo/yr. 



(5) 



(iv) Old wide binaries. Another limit comes from the 
existence of old wide binaries, since for a very mas- 
sive galactic nucleus the galactic orbits would have been 
much smaller than at present. Therefore the density of 
stars would have been much larger and the dissolution 
time of binaries due to stellar encounters correspondingly 
shorter. From a list of 11 well observed old wide binaries 
one finds a limit on steady mass loss @ 



M < O(lO)M0/yr. 



(6) 



In conclusion, we have four different methods which all 
give a bound on — M between 0(1) and 0(10) Af©, so we 
write 



Af<e Mo/yr, 



(7) 



and we expect e '^ 1, unless one cannot find a way, per- 
haps with some fine tuning, to relax the most stringent 
bound I0J). Anyway, given that we have three more in- 
dependent limits on e, we see that we cannot stretch 
the value of e beyond, say, e ~ 10. Furthermore, inde- 
pendently of the technical details leading to the above 
bounds, it is easy to understand qualitatively why a 
bound on — M of order IMo/yr emerges. The total mass 
of the galactic disk plus bulge and spheroid is estimated 
to be 9 X 10^°Mq, while a lower bound on the age of the 
Galaxy is provided by the age of its oldest globular clus- 
ters, which is 1.2 x lO"'^" yr. Therefore a mass loss rate 
due to GWs of the order of a few solar masses per year 
implies a total mass loss, over the age of the Galaxy, com- 
parable to its present mass. It is not surprising that, at 



this level, one finds important dynamical effects related 
mass loss. If one streched e up to values of order 10, one 
should admit that the mass that the galactic disk has 
lost to GWs during its history is larger than its present 
mass, so over 50% of the original mass would have been 
lost to GWs. 

We also remark that these are just upper bounds and, 
as far as the above arguments are concerned, there is no 
physical reason that suggests that they might be almost 
saturated, so they should not be taken as an indication 
of a plausible value of M. 



III. IMPLICATIONS FOR GALACTIC BURSTS 
OF GWS 

We now compare this bound, which is generic and 
holds whatever is the physical origin of the mass loss 
of the Galaxy, with the energy carried away by GWs, if 
there are GW bursts with a typical rss amplitude h^ss 
and a rate TZ. Following Ref . jj] , we consider a Gaussian 
waveform of duration r, given by 
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whose Fourier transform is 

Mf) = /i.ss(2^r2)i/4 exp(_7rV2/') • (9) 

(The same analysis can be repeated with sine-gaussian 
waveforms, with similar results). We first consider a wave 
coming from optimal direction with h^[t) = h[t) and 
/ix(0 = 0- Fi'om the standard expression of the energy 
flux. 
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where dA — r^dil and r is the distance to the source, we 
get that the total energy radiated by such a burst, 



E = 47rr' 



dE r^c^hl. 



df AGt^ 



(11) 



The average over the square of the pattern function of the 
interferometer gives the usual factor 2/5 which, following 
the conventions of the LIGO collaboration, we include in 
the definition of /irss- Then, the energy carried by a burst, 
averaged over the arrival direction and the polarization, 
is 



{E) = 



r-^c^K' 



AGt^ 
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8 kpc 



10-19 Hz 
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1 ms 



(12) 



where in the second line we normalized r to the distance 
to the galactic center. We see from the above that a burst 




FIG. 1: The limit from galactic dinamics with e = 0.5 (solid 
line) and e = 10 (dashed), compared to the bound from the 
LIGO S2 run, when r = 1 ms. LIGO data taken from Fig. 12 
of Ref. Q. 



of duration 1 ms, with h^-ss — 10^^^ Hz^ ' , carries away 
about 3 X 10^^ solar masses in GWs, if it comes from 
a source located at typical galactic distances. Of course 
shorter bursts are more energetic, since they extend in 
frequency space up to /max ~ I/t. In an astrophysi- 
cal context it makes sense to consider bursts as short as 
r ~ 0.1 ms, corresponding, for the gaussian wavepacket 
® ; to /max of Order of a few kHz, which indeed are the 
shortest burst searched for in Refs. 0,l2|- If 7?. is the rate 
of these bursts, the associated rate of mass loss of the 
Galaxy is M = -7^ x {E)/c^, with {E) given by eq. ((1^ . 
Combining this with eq. {TJ we get 



7^< 8.0e- 



events fSkpc^ flQ-^^Rz"^^^^ 
day 



1 ms 



.(13) 
where now r is an average distance scale characterizing 
the population of sources [la| . In Fig. ^ we show this 
bound, in the [TZ, /irss) plane, for r = 1 ms, r = 8 kpc, 
and we compare it to the experimental bounds obtained 
by the S2 LIGO run, for the same value of r. The solid 
line corresponds to e = 0.5, and the dashed line to e = 10. 
In Fig.|21we perform the same comparison for t = 0.1 ms. 
We see that the bound from galactic dynamics is a 
significant one, especially for small r, and indeed for 
T = 0.1 ms is quite stronger than the LIGO S2 re- 
sult. With the on- going and future high-sensitivity runs 
by LIGO and VIRGO one can expect that, even for 
T = 0.1 ms, one will go beyond the small-Zij-ss portion 



of the limiting curve determined from galactic dynamics, 
while at larger h^-^g the limit from galactic dynamics will 
remain the dominant one. If a statistically significant 
rate in excess of the value obtained from eq. 113|) with 
e ~ 1 and r ~ 8 kpc were found and if, from the energetic 
of the events, one concluded that they originated in our 
Galaxy, then one should study the possibility of relax- 
ing somehow the bound on e (which anyway should be 
possible at most by one order of magnitude) . Otherwise, 
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FIG. 2: The same as Fig. 1, with r = 0.1 ms. 



one could consider the possibility that, rather than hav- 
ing a homogeneous distribution of sources in the galactic 
disk, the signal could be due to a single (or a few) source 
at a close distance from us, which emits repeatedly GW 
bursts [lOl lli^|. Therefore, eq. H13(l can give useful clues 
as to the spatial distribution and possibly the origin of 
the sources. 

Conversely, if GW detectors should find that the bound 
that we have discussed is saturated, this would mean that 
the effect of GW emission on some astronomical observ- 
ables, such as the radial motion of the LSR, is just of 
order 10-20 % of the present observational uncertainties. 
This would imply that, with an increase in the accu- 
racy of astronomical observations, as is expected with 
the GAIA mission (see e.g. Ref. Q), one could be able 
to single out the effect of GW emission on the dynamics 
of the Galaxy, providing a form of indirect detection of 
GWs. 
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